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Abstract 

With the advent of the recent measurements in neutrino physics, we investigate the 
role of high-energy neutrino flux ratios at neutrino telescopes for the possibility of de- 
termining the leptonic CP-violating phase 5 and the underlying pattern of the leptonic 
mixing matrix. We find that, even considering the current uncertainties on the neutrino 
mixing parameters, the flux ratios show a dependence of O(10 %) on the CP-violating 
phase, and for optimistic uncertainties on the flux ratios less than 10 %, they can be 
used to distinguish between CP-conserving and CP-violating values of the phase at la 
in a non- vanishing interval around the maximal value |<5| = ir/2. We also show that 
the correlations existing among the flux ratios can be useful to discriminate different 
classes of models based on tri-bi-maximal and bi-maximal mixing patterns. 
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1 Introduction 



Ever since 1998 with the first results of the Super-Kamiokande experiment on atmospheric 
neutrinos pQ, there is strong evidence for neutrino oscillations as the main mechanism for 
neutrino flavor transitions. Indeed, the fundamental neutrino parameters such as the two 
neutrino mass-squared differences Am| x and Am^ as well as the three leptonic mixing angles 
6*23; 6*12; and #13 are now determined with better and better accuracy. Recently, the third and 
last mixing angle #13 has been measured by the Daya Bay experiment [2] , but also the Double 
Chooz, MINOS, RENO, and T2K experiments have made important contributions [3H6]. All 
five experiments indicate that the value of the third mixing angle is around nine degrees. 
Thus, the value of this mixing angle is relatively large. This means that the tri-bi-maximal 
(TBM) mixing pattern for leptonic flavor mixing has essentially been ruled out as a zeroth- 
order approximation, and similarly for the bi-maximal (BM) pattern. Therefore, appropriate 
corrections must be taken into account to reconcile them with the experimental data. To this 
end, there remain two quantities for massive and mixed neutrinos that can be determined 
with neutrino oscillations: the sign of the large mass-squared difference sgn(Am| 1 ) and the 
Dirac CP-violating phase 5. The latter always appears in combination with #13 in the leptonic 
mixing matrix as sin($i3) exp(±i5) and a non-zero value of #13 basically means that it is 
possible to determine 5. The best options to measure S are provided by accelerator neutrino 
experiments (e.g. NOvA and NuMI), future superbeams, beta beams, or even better a so- 
called neutrino factory [7]. However, it will take a long time before such experiments are 
realized. In this work, we will therefore consider an alternative approach to determine § by 
measuring neutrino flux ratios at neutrino telescopes. The generic example of such a neutrino 
telescope is the IceCube experiment [Hj, which exists and (although not straightforwardly) 
has the potential to measure neutrino flux ratios [9] . The dependence of neutrino flux ratios 
on the neutrino mixing parameters has been previously addressed in the literature where 
the focus has been on the dependence on 6*23 and 6*13 [T0TfT2"] . on 5 [13], and on new physics 
effects pp. 

The aim of this work is to illustrate the potential of neutrino telescopes to detect the 
CP-violating phase 5 through the measurement of neutrino flux ratios after the first mea- 
surements of the mixing angle #13. After introducing the general definition of flux ratios 
in Section [2l we explicitly show the 5 dependence of these flux ratios and discuss the role 
of the updated mixing parameter uncertainties in Section [31 This will help in identifying 
the main source of error in the theoretical estimate of the flux ratios, assessing especially 
the role of #13. Assuming that the experimental flux ratios will be measured within a given 
precision, we estimate the capability of the neutrino telescopes to distinguish, at 2a, between 
CP-conserving and CP- violating values of 5. Then, in Section HI we study for the first time 
to what extent flavor models based on TBM and BM mixing can be differentiated at neu- 
trino telescopes. Finally, in Section we summarize the results and present our conclusions. 
For our numerical evaluations, we adopt the best-fit values and uncertainties for the mixing 
angles in the normal neutrino mass hierarchy (NH) quoted in Ref. [TO], which are presented 
in Tab. [TJ for completeness. Note that, although the neutrino flux ratios are not sensitive 
to the neutrino mass-squared differences, a somewhat dissimilar behavior between NH and 
inverted hierarchy (IH) could be drawn if the allowed ranges for the parameters were sizeably 
different. However, this is not the case for the 3a ranges considered in this analysis. 
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Parameter 


Best-fit value 


la range 


2cr range 


3cr range 


sin 2 Ou/10- 1 


3.2 


3.03 - 3.35 


2.9 - 3.5 


2.7-3.7 


sin 2 e 13 /io-' 2 


2.6 


2.2 - 2.9 


1.9 - 3.3 


1.5 - 3.6 


sin 2 023/1O- 1 


4.9 


4.4 - 5.7 


4.1 - 6.2 


3.9 - 6.4 


6/ir 


0.83 


0.19 - 1.37 


[0,2tt] 


[0,2tt] 



Table 1: Results of a global three neutrino oscillation analysis f76] / in terms of best-fit values, 
la, 2a ; and 3cr ranges for NH only. See also Ref. J7?| / for another global analysis. 



2 Neutrino flux ratios 

The averaged neutrino oscillations probabilities can be written as (see e.g. Ref. [TO] ) 
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i=i 



where the quantities U a i are elements of the leptonic mixing matrix U. Starting from the 
flux ratios at a source given by (fP Ue : </>° : <f)® T , the neutrino fluxes arriving at a detector are 
sensitive to neutrino oscillations in vacuum and are then computed as 

k= E « p ^>= E <Ei^iW- ( 2 ) 

For example, in the case : : 0^ = 1:2:0, the neutrino fluxes are stemming from 
pion-beam sources (ttS), whereas in the case : 0° : = 0:1:0, they come from 
muon-damped sources (fiDS). Then, the three-flavor neutrino flux ratios are defined as 
follows: 

Rafl = ^- 0) 
<Pi/p 

Note that the three neutrino flux ratios R elx , R eT , and R^ T are not independent of each other, 
since two of them will give the third one. For example, using i? eAt (i? er ) _1 i? MT = 1, we showed 
in Ref. [10] that, up to second order in small quantities (to be defined below), the following 
sum-rule holds: R efJL — R er + R^ ~ 1, which enable us to eliminate one ratio in favor of 
the other two. Given also the simplicity of their analytical expressions, we will choose to 
examine first the relevant features of two neutrino flux ratios R eii and R^. Notice that a 
more accessible variable from the experimental point of view is given by 

R = JZ+^Z = R,ie T+R7e = (i ^ rl l+(i£F ' (4) 

which is obviously related to the "fundamental" R a p- We will take into account R in Sec- 
tion I3.3[ when dealing with the issue of achievable precision for the flux ratios at neutrino 
telescopes. It is important to stress that the various flux ratios are experimentally accessible 
at different neutrino energies. For example, IceCube has an energy threshold of 100 GeV 
for detecting muon tracks, and ~ 1 TeV for detecting electron- and tau-related showers. At 
higher energies, above ~ 1 PeV, electron-related electromagnetic showers and the tau-related 
hadronic showers can also be distinguished. Here, we do not consider such details and assume 
to work in the appropriate energy regime, where all relevant R a p can be measured. 
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3 Investigation of neutrino flux ratios 



3.1 Estimate of the uncertainties for pion-beam sources 

Since the exact formulas for the neutrino flux ratios are quite cumbersome, we prefer to 
present our results expanding these ratios in the small parameters #13 (although there is a 
best-fit value for this parameter, it is still small), #23 = #23 — 7r/4, and 612 = 612 — #12, #12 
being the best-fit value for 9\2. No restrictions have been applied to the CP- violating phase 
5, which means that the following formulas are valid to all orders in 5. In our discussion, it 
is enough to consider the expansions up to first order for R efl and R eT and second order for 
R^ T (see Ref. [10] for a detailed discussion on first- and second-order approximations from 
maximal mixing), we obtain 

1 + 3 cos(5) sin ( 4 £ 12 ) 13 _ I s i n 2 (20 12 ) 5 23 , 

(5) 

1 + 2 cos 2 (<5) sin 2 (2# 12 ) 9 2 13 + 2 cos(*) sin(4# 12 ) 9 13 5 23 + [cos(4# 12 ) + 7] * 2 3 . 

At this order of perturbation theory, R eT = R e ^. Notice also that no contributions from 
612 appear at the perturbative level considered above, so we expect that, given the current 
best-fit value on the mixing angle 612, the impact of 5 12 is negligible. Inserting the best-fit 
values from Tab. [JJ into Eq. (JSJ), we can evaluate the numerical weight of each term in the 
expansions: 

~ l + 0.5cos((J) 0i3- 1.3*23 + O(<5?0» 

(6) 

R^ T ~ 1 + 1.7 COS 2 (5)0 2 3 + 1.3 COs(*) ^13*23 + 6.3 * 2 3 . 

Numerically, *i 2 ~ 0.05, #13 ~ 0.2, and *23 ~ 0.15. Since *23 always appears with the largest 
coefficient, we can easily deduce that the main uncertainty for the flux ratios comes from 
the current error on the mixing angle ^23- The contribution of 6* 13 is modulated by *, so we 
expect no impact from it for cos(<5) ~ 0. 

We verified all these statements by computing the R a pS as functions of *, using exact 
expressions for the neutrino flux ratios. For the sake of illustration, we show in Fig. [1] the 
results obtained for both R e ^ and R^ T . In the upper plots, we show the behavior of the 
two flux ratios when 13 is varied inside its 3a range, whereas in the lower plots, we show 
the same for 823, with the unshown parameters fixed to their best-fit values given in Tab. [TJ 
These plots confirm the conclusions drawn from Eq. (JSJ). Although not explicitly shown, the 
impact of 6*12 is completely negligible for R^r, whereas for R etl , it is numerically comparable 
to that of 013. 

As a next step, we evaluate all possible allowed values of the neutrino flux ratios con- 
sidering the simultaneous variations of all parameters within their respective uncertainties 
as specified in Tab. [TJ The results are shown in Fig. [2] (upper plots), where we also show 
the present correlation between the two neutrino flux ratios (lower plot). In these plots, the 
bands indicate allowed regions where all parameters are varied within their respective 3a 
ranges, whereas the solid and dashed curves have been obtained from the 2a and la ranges, 
respectively. Figure |2] deserves some more comments. Let us first consider the allowed re- 
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Figure 1: Variations of R e ^ (left plots) and R^ T (right plots) as functions of 5 when only the 
uncertainty on #13 is taken into account (upper plots) or that on 823 (lower plots). 



gions within the dashed curves, which can be used to illustrate a possible future situation 
when all mixing parameters (especially 6*23) will be known with better accuracy than today. 
In this case, the dependence of the neutrino flux ratios on 5 is quite pronounced, since i? e/i 
has maximum and minimum values corresponding to the CP-conserving values 5 = 0, ±7r, 
respectively; at 5 = ±7r, R efl G [0.83,0.97] and at 5 = 0, R efl G [0.96,1.13]. For R^ T , the 
5 dependence is a bit more complicated, so that the flavor ratio assumes similar values in 
[1.03, 1.09] for 5 = 0, ±vr, but R^ G [1.00, 1.03] for 5 = =br/2. 

If the errors on the mixing angles were larger, the dependence on 5 would be less pro- 
nounced, but still enough to be used to extract information on it (see below). Globally, we 
obtained the following ranges for the flux ratios: 

R eii G [0.85,1.18], iV G [1.00,1.21]. 

The correlation between the two neutrino flux ratios shown in the lower plot of Fig. [2] allows 
one to make some consistency checks. In fact, working in the absence of exotic processes 
(like neutrino decays), two simultaneous measurements of R etli and R^ T should fall in these 
regions. For example, R efl ~ 0.8 must be accompanied by R^ T ~ 1.1 to give \5\ ~ ir. Indeed, 
it is interesting to note that the authors of Ref. [17] found a possible hint in favor of 5 ~ ir. 
Thus, given all other best-fit values, these results - if they are confirmed - would imply 
{Ren, Rfj, T ) = (0.93, 1.04), which is indicated by a "star" in the lower plot of Fig. [2l 
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Figure 2: Total uncertainties on R efl (upper-left plot) and R^ (upper-right plot) as functions 
of 5 when all other mixing parameters are varied within their 3a errors. In the lower plot, 
we show the correlation between R efl and R^. Bands indicate the allowed regions when all 
parameters are varied in their 3a ranges, whereas the solid and dashed curves are obtained 
from the 2a and la ranges, respectively. The "star" represents the flux ratios at the best-fit 
point of the neutrino parameters. Note that 5 is left free in the interval [— ir, 7r]. 



3.2 Discussion on muon-damped sources 



A similar study can be carried out for a muon-damped source, for which A 



'v e ■ ■ 4>l T - 

0:1:0. Up to first order in the small parameters d~ij, the relevant neutrino flux ratios read 



R 



4sin 2 (20 12 ) 64sin(40 12 ) 4 cos(5)[-9 + cos(40 12 )] sin(40 12 ) 

+ 7Z~. — mo 012 7Z ; , An m, #13 



7 + cos(40 12 ) [7 + cos(40 12 )p 
8[-9 + cos(40 12 )]sin 2 (20 12 ) . 

H — — „n 02Z 



[7 + cos(40 12 )] 5 



R 



[7 + cos(40 12 )] 2 
4 cos(5) sin(46'i 2 ) 



7 + cos(40 12 ) 13 ' 7 + cos(40 12 
Inserting the best-fit values of the mixing parameters, we obtain 

Rep = 0.6 + 1.1 5i2 + 0.7 C0S(5) 013 -1.7*23 + 0(<SJ 



8 sin 2 (20 12 ) 



5 



23 



R 



1.0 - 0.4 cos(5) 0i3 + 0.9 S 23 + 0{b\ 



(7) 



(8) 
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where the uncertainty from 5 12 is negligible at this order in R^. Some important differences 
arise when comparing the formulas in Eq. (jSJ) with the corresponding ones in Eq. (JSJ). In fact, 
R efJL now obtains a leading dependence on 5 12 , meaning that this flux ratio is more sensitive 
to the uncertainty of the mixing angle 612 than the corresponding ratio for the 1:2:0 case. 
At the same time, R^ T is corrected by the standard unit value by linear terms in #13 and #23, 
which means that the fluctuations due to the current uncertainties on the mixing angles are 
much more pronounced than before. These considerations are based on simple analytical 
formulas that have been verified numerically, computing the exact flux ratios as functions 
of 8. As for the previous case, we found that the uncertainty on #23 is the dominant source 
of error, and for R efl , the error induced by 6*12 is a bit larger than the one given by 6*13. 
Besides these differences, the qualitative behavior of the neutrino flux ratios as functions 
of 5 is similar to the one shown in the upper plots of Fig. |2l so we have not included the 
corresponding figures. On the other hand, the correlation between the two flux ratios is 
quite different from the previous case, and we display it in Fig. [3j As can be observed, there 
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Figure 3: Correlation between R ejl and R^ T for the 0:1:0 case obtained when the mixing 
parameters vary within their la, 2a, and 3a ranges. Note that 5 is left free in the interval 

[-7T,7r]. 

exists a "negative" correlation, extending for R e>1 G [0.27,0.92] and R^ T G [0.92, 1.42]. 

3.3 Potential to measure a non-vanishing CP-violating phase 

It is clear that all previous analytical considerations should be supported by experimental 
feasibility in measuring the neutrino flux ratios. Following Ref . [H] , we study the potential 
of neutrino telescopes to distinguish between CP-conserving and CP-violating values of 8, 
using the possible measurements of R and R eT and also the capability to differentiate between 
nS and [iDS sources. This is a difficult task, as was outlined for example in Ref. [9], 
especially if the uncertainties on the considered flux ratios are relatively large. For the 
sake of completeness, we present the explicit dependence on 5 for the flux ratios R and R eT 
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computed at the best-fit values of 0y for both nS and f^DS sources: 



R 



R »DS 



1.99 + 0.03 cos(<5) 
1 

0.61 + 0.03 cos(cT) 



RI mDS 



19.1 + cos(<J) 



18.7 - 0.47cos(5) - 0.21 cos(2<5) ' 

8.94 + cos(5) 
" 14.3 + 0.02 cos(cT) - 0.21 cos(25) ' 



(9) 



The coefficients of the cos(<5) terms are quite small and relevant for the CP-conserving 
values. In addition, there are no preferred sources to be used in the determination of the 
CP-violating phase, since the 5 dependence in the flux ratios built from nS and fiDS is not 
really dissimilar. 

To evaluate how 5 = can be distinguished from other non-vanishing values, we build a 
simple x 2 function 



X 



i? ex P-i?(%,5) 



or 



Rl* p — R eT (9ij, S) 



(10) 



where R cxp and R! 1 ^ are the ratios R and R er evaluated at the best-fit points for 9ij as given 
in Tab. [U We consider the possibility of having the sum of both contributions from irS and 
fiDS sources. In the minimization procedure, we marginalize over all parameters but 5. As 
an illustrative (and favorable) example, we show in Fig. H] the case where S = ir/2; the three 
lines refer to nS sources (dashed curve), fiDS sources (thin solid curve), and the sum of 
both (thick solid curve). In the \ 2 function, we assume an optimistic 5 % error on the flux 
ratios. As expected, we see that the combination of both irS and fiDS sources gives the 
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Figure 4: \ 2 ~ Xmin as a function of 5, after minimizing over all non-displayed parameters. 
The three curves refer to nS (dashed curve), fiDS (thin solid curve), and the sum of both 
(thick solid curve). 

best resolution for S. In particular, the x 2 function does not touch the CP-conserving values 
5 = (0, 7r) at 2a, so maximal CP violation can be distinguished from CP conservation at a 
good confidence level. This exercise can be repeated for every input value of 5; at a given 
confidence level and flux ratio uncertainty, there exist a range of phases for which the \ 2 
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function does not touch 5 = (0, 7r); the fraction of such points over the whole [0, 7r] range 
defines a CP fraction that serves to illustrate the goodness of the neutrino telescope to access 
CP violation. The result of such an exercise is reported in Fig. [51 where we display the 2a 
CP fraction in terms of the uncertainty <jr = <jR eT = a, ranging from (the unrealistic) 1 % 
to (the more realistic) 10 %, obtained using the \ 2 function defined in Eq. ( fTOj) . We have 
considered four cases, which correspond to ttS + fiDS (solid curve), fiDS (dashed curve), 
and 7i S (dot-dashed curve); the long-dashed curve refers to the case in which all mixing 
parameters but 5 are fixed to their best-fit values and all neutrino sources are taken into 
account. It is clear that the best performance of a neutrino telescope is reached when the 

0.7 
0.6 

o 

3 0.4 
* 0.3 
U 0.2 
0.1 



"0".00 0.02 0.04 0.06 0.08 0.10 

cr (%) 

Figure 5: 2a CP fraction estimate as a function of the common flux ratio uncertainty. The 
four cases correspond to nS + fiDS (solid and long-dashed curves), fiDS (dashed curve), 
and irS (dot-dashed curve). 

maximum amount of information is collected (i.e., the sum of the sources). However, such 
performances are limited to a < 0.065 and do not appear to be particularly encouraging. 
Even considering infinite precision on the mixing angles, the CP fraction is non-vanishing 
only for a < 0.07. This is the consequence of the mild (^-dependence on the flux ratios in 
Eq. ©. 

4 Checking flavor models 

The correlation discussed in Section [3] can also be useful, at least to some extent, to 
discriminate among different classes of flavor models. It is clear that, whatever the leading- 
order approximation of the neutrino mass matrix is, every successful model should overlap to 
some portion in the (-R e/i , i? AtT )-plane, since they have to reproduce the experimental values 
of the various mixing angles. However, it could happen that a specific model predicts a 
well-defined correlation between the two neutrino flux ratios that can be checked at neutrino 
telescopes to favor or disfavor the presumed neutrino mass texture. 

To illustrate this point, we consider the prediction for the mixing angles of a typical Aa 
model for the TBM pattern pHEH] and a typical S 4 model for the BM pattern [2OH220 

both TBM and BM, sin 2 (6> 2 3) = 1/2, 9 13 = 0, and the mixing angle 6 12 is given by sin 2 (6>i 2 ) = 1/3 
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Following Ref. [23], the results of the next-to-leading-order corrections to the TBM pattern 
produce 



(sin 2 # 23 )tbm = 2 + Re ( c 2s) £ + ~j= 



Re(c^) - v / 2Re(c y 



23) 



(sin 2 12 ) TBM = ^ - ?Re(c^ 2 + cf 3 ) £ + ^ Re(c^) £ , W 



fsin 6 1 



3 3 
1 



13JTBM 



6 



3V2 (c e 12 - c e 13 ) + 2V3 (V2 Ci 3 + c: 



where the coefficients are random complex numbers with absolute values of order 1 and 
the ij indices are the corresponding entries of unitary matrices diagonalizing the charged 
lepton and neutrino mass matrices at the next-to-leading order. The small parameter £ has 
been chosen in such a way to maximize the success rate to reproduce all three mixing angles 
inside their corresponding 3a ranges, and it turns out that this happens for £ ~ 0.07. A 
similar approach to the corrections for the BM pattern gives 

111 1 
(sin 2 9 23 )bm = g > ( sin2 Mbm = g - T7| Re ( c i2 + c is) f > (sin# 13 ) B M = |cf 2 - c^l £ , 

(12) 

where all coefficients have the same meaning as before, but instead £ ~ 0.163. Starting from 
the previous expressions, we compute the values of the mixing angles generating (9(10 5 ) 
with a Gaussian distribution around 1 with variance 0.5. The results are summarized 

I LJ 

in Fig. [6j In the left plots, we show the result of our parameter scan for the TBM pattern. 
As can be observed, the next-to-leading order corrections are enough to generate neutrino 
flux ratios in almost the whole plane. The points tend to cluster close to the lower limit of 
R^r as a consequence of the fact that larger values can only be obtained by huge deviations 
for #23 from being maximal and \8\ close to ir (see the upper-right plot in Fig. |2j). For the 
BM pattern, i.e., for the right plots, the situation is quite different, essentially because 623 
is kept fixed to maximal mixing (but we checked that introducing 0{^ 2 ) corrections do not 
affect our results too much). We observe that the strong correlation between R efl and 
can be used to rule out the BM pattern as a leading-order approximation if the neutrino flux 
ratios sufficiently deviate from 1. 

To understand the range of values of R efJi and R^ T for these classes of models, one has to 
expand the probabilities around the leading-order result (TBM or BM) and take into account 
the signs and order of magnitudes of the shifts. This has been discussed in the Appendix to 
which we refer the interested reader. 



5 Summary and conclusions 

In this paper, we have carefully analyzed the potential of neutrino telescopes to access 
the leptonic CP-violating phase 5 and their capability to distinguish among different classes 
of flavor model for leptonic mixing. We have derived analytical perturbative expansions for 
the flux ratios R a p up to first (and second) order in small parameters, explicitly showing 
their dependence on 5, for both kinds of high-energy neutrino sources, i.e., irS and fiDS 
sources. It turns out that the uncertainty on 823 affects the global (theoretical) error on the 

and siii 2 (#i 2 ) = 1/2, respectively. 
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Figure 6: Results of parameter scans for the TBM (left plots) and BM (right plots) patterns 
in the (R efl , R^ T ) -plane. 

flux ratios the most. Considering at the same time both kinds of sources, we have shown that 
a 10 % of CP fraction can still be obtained with a 5.5 % uncertainty on R a p. The correlation 
existing among various flux ratios can be used to discriminate, to some extent, models of 
neutrino masses based on discrete symmetries. In particular, we have found that textures 
of bi-maximal type at leading order (with vanishing next-to-leading order corrections on 
#23) can be ruled out if R^ is larger than 1.2 and R efl is sufficiently different from the 
interval [0.4,0.6]. Finally, we urge the IceCube collaboration to measure the neutrino flux 
ratios, since such a measurement could provide the first hints on the value of the leptonic 
CP-violating phase. 
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A Expansions of flux ratios 



A.l Expansion of the flux ratios around TBM 

In this appendix, we expand the R efl and R^ T flux ratios around the TBM ansat 
for sin 2 (#23) = 1/2, #13 = 0, and sin 2 (#i 2 ) = 1/3. For pion-beam sources, we find that 



(#L)tbm 



1 + ~Y C0S ( 5 ) 013 - ^23 , 



(13) 
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I COS (5) <? 2 3 + \/2cos(<5) e 13 5 2 3 + 75 



Similarly, for muon-damped sources, we obtain 



4 3fi 22 88 

« 5 )tbm = - + -v / 2 5 12 + -v / 2cos(5)^ 13 --5 23 



(14) 



Let us consider the results of Fig. [6] for TBM-ttS 1 (upper- left plot). From the expression 
of (-R^J)tbm, we see that the maximal value is obtained when 623 < and cos(<5) = 1, 
whereas the minimum is reached for opposite values (#13 is always larger than zero). Given 
the maximal allowed values for these variables, we obtain (-R^)tbm £ [0.7, 1.2], as shown in 



the scatter plot. Similar arguments can be invoked to explain the range for (R 
for the flux ratios from fiDS. 



irS-) 

[IT 



/TBM 



and 



A. 2 Expansion of the flux ratios around BM 

Here, we expand the R e ^ and R^ T flux ratios around the BM ansatz, i.e., for sin 2 (6* 12 ) = 
sin 2 (#23) = 1/2 and 13 = 0. For pion-beam sources, we find that 

(Re^)BU = 1 — 7^ ^23 , 

(15) 

W)bm = l + 2cos 2 (5)fl 2 3 + 6<5 2 3 . 

It is interesting to observe that (R^)bm ~ (^«)bmi as suggested by the upper-right plot 
in Fig. EJ Similarly, for muon-damped sources, we obtain 



(-^e/f 5 )bm - 3 - y 5 23 



« S )BM = 1 + ^23 



(16) 

3 
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